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Abstract 

In this article we consider the problem of giving a robust, model- independent, lower 

bound on the price of a forward starting straddle with payoff \Fti — Ft,, \ where < To < Ti . 

CO Rather than assuming a model for the underlying forward price {Ft)t>o, we assume that call 

prices for maturities To < Ti are given and hence that the marginal laws of the underlying 

are known. The primal problem is to find the model which is consistent with the observed 

, call prices, and for which the price of the forward starting straddle is minimised. The dual 

f-) . problem is to find the cheapest semi-static subhedge. 

.^ Under an assumption on the supports of the marginal laws, but no assumption that 

the laws are atom-free or in any other way regular, we derive explicit expressions for the 
coupling which minimises the price of the option, and the form of the semi-static subhedge. 
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We consider the problem of constructing martingales with two given marginals which minimise 
the expected value of the modulus of the increment. This problem has a direct correspondence 
to the problem in mathematical finance of giving a no-arbitrage lower bound on the price of 
an at-the-money forward starting straddle, given today's vanilla call prices at the two relevant 
maturities. There is also a related dual problem, which is to construct the most expensive 
semi-static hedging strategy which sub-replicates the payoff of the forward starting straddle for 
'^ any realised path of the underlying forward price process. Under a certain assumption on the 

fvj distributions we solve the primal and dual problem and demonstrate that there is no duality 

^- gap. 

^^ The results of this article complement previous results by Hobson and Neuberger j20j. In 

cn that article, the authors solved the analogous problem of constructing no-arbitrage upper price 

. . bounds and semi-static super-replicating hedging strategies for the forward starting straddle. 

_^ Hobson and Neuberger [20j give some examples, but the main contribution is an existence result, 

S^ and a proof that the primal and dual problems yield equal values. 

Returning to the lower bound case, it follows from results of Beiglbock et al [2| that there 
exists a solution to the primal problem and that there is no duality gap. However, Beiglbock et 
al [2j give an example to show that the dual supremum may not be attained. Our contribution 
is to show that, under a critical but natural condition on the starting and terminal law, there is 
an explicit construction of all the quantities of interest. (A parallel construction gives an explicit 
form for the optimal martingale in the upper bound setting, but in that case the condition on 
the measures is less natural.) 
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There has been a recent resurgence of interest in problems of this type, in part motivated by 
the connection with finance, see the survey article by Hobson [T7], and in part motivated by the 
connections with the optimal transport problem, see [2l[3l[T5]. A first strand of the mathematical 
finance literature is concerned with robust pricing of particular exotic derivatives, for example 
lookback options [Ml [Hj , barrier options [3 [H [9] , Asian options [12] , basket options [HI [10] 
and volatility swaps [3 HU [TH] . This strand of the literature often makes use of connections 
with the Skorokhod embedding problem, but in some senses the problem here is simpler in that 
the option payoff only depends on the joint law of {Ft^, Fti) and is otherwise path- independent. 
For this reason the full machinery of the Skorokhod embedding problem is not needed, although 
it can still help with the intuition. A second strand of the mathematical finance literature on 
robust pricing [HJ HU JU [l3] considers consistency between options and no-arbitrage conditions. 

The optimal transport literature (see Villani |25j ) is concerned with the cheapest way to 
transport 'sand' distributed according to a source measure fj, to a destination measure i^. Such 
problems are also motivated by economic questions and the transport of mass becomes an issue 
of the optimal allocation of economic goods. Unsurprisingly, questions of this type were of par- 
ticular importance to mathematicians working on questions of efficiency in planned economies, 
see for instance Kantorovich [^ . The Kantorovich relaxation of Monge's original problem goes 
back to his use of linear programming methods and as such, the Lagrangian approach taken 
in this paper is most natural. With respect to the classical optimal transport literature, the 
novelty of the current problem, as elucidated in Beiglbock et al [2] and developed in Beiglbock 
and Juillet [3] and Henry-Labordere and Touzi |15| is to add a martingale requirement to the 
transport plan, which is motivated by the idea that no-arbitrage considerations equate to a con- 
dition that forward prices are martingales under a pricing measure. Both [3j and [15] consider 
the problem of minimising the martingale transport cost for a class of cost functionals , but 
the payoff |y — x| is not a member of this class. In contrast, here we focus exclusively on the 
payoff |y — 3;|. This functional encapsulates the payoff of a forward starting straddle, which is 
an important and simple financial product which fits into the general framework, and it is the 
original Monge cost function in the classical set-up. For these two reasons this cost functional 
is of significant interest. 

2 Motivation and preliminaries 

Let X, Y be real-valued random variables and suppose X ~ ^u and Y ^ u. We say that the 
bivariate law p is a martingale coupling (equivalently a martingale transference plan) of X and 
Y, and write p G M.{fJ,, v), if p has marginals p, and u and is such that J (y — x)p{dx,dy) = for 
each X. For a univariate measure C define C(^ via C^(x) = J{y — x)^({dy). It is well known (see, 
for instance, Strassen[24J) that Ai{p, v) is non-empty if and only if p is less than or equal to v 
in convex order, or equivalently /x and v have equal means and C^{x) < C,^{x) for all x G M. 
Then, under the assumption that M.{p, v) is non-empty, our goal is to find 

r{p,u):= inf E[\Y-X\]= inf [ \y - x\ p{dx , dy) . (2.1) 

The financial significance of this result is as follows. Let F = {Ft)t>o denote the forward 
price process of a financial asset. Let Tq and Ti be two future times with < Tq < Ti. A well 
known argument due to Breeden and Litzenberger |4J, shows that knowledge of a continuum 
of call prices for a fixed expiry is equivalent to knowledge of the marginal law at that time. 
Suppose then that a continuum in strike of call prices are available from a financial market and 
hence that it is possible to infer the marginal laws of F at time Tq and time Ti. Suppose these 



laws are given by /i and v. The problem is to minimise E[|i*Ti ~-^Tol] over all martingale models 



for F with the given marginals, i.e. to find (2.1). We will call this problem the primal problem. 



Conversely, suppose we can construct a trio of functions (V'O;^!) f^) such that 

\y - x\ > i^iiy) - M^) + 5{x){x - y) Vx, y G M. (2.2) 



Then E[\Y - X\] > E[V'i(^)] - E[V'o(^)] = j ^i{y)y{dy) - fjMx)fi{dx). It follows from (2.2) 



that if S is the set of trios of functions (^O) V'Ij ^) such that (2.2 ) holds then 



E[\Y - X\] > sup <j / MvMdy) - / M^)Kdx) } =: V{fi, v). (2.3) 

(Vo,V'i,<5)e5 



We will call the problem of finding the right-hand-side of (2.3) the dual problem. 

Again there is a direct financial interpretation of ( |2.3| ) . Under our assumption that a contin- 
uum of calls is traded, the European contingent claims V'i(-Pti) and '0o(-^To) can be replicated 
with portfolios of call options bought and sold at time zero. Moreover, if the agent sells forward 
5{Fxq) units over [To,Ti] then the gains are 5{Fj'q){F'j'q — F^^)- Combining the two-elements 
of the semi-static strategy ([TTl Section 2.6]) consisting of calls and a simple forward position 
yields 

MFt,) - MFto) + S{Ft,){Ft, - Ft,) 



which corresponds to the right-hand-side of (2.2). If (2.2) holds then the semi-static strategy is 
a subhedge for the payoff | Ft, — Ftq \ ■ 



It is clear that if (2.2) holds then we must have ipi{x) < ^po{x) and that if we want to 



find ipi to maximise the right-hand-side of (2.3) then we want ^o as small as possible. As 



a result a natural candidate for optimality is to take ^o = V'l and the problem of finding a 
trio {^jJQ,^pl,6) reduces to finding a pair (V'j^)- We let S be the set of pairs {ip,5) such that 
|y — 2;| > ip{y) — ip{x) + S{x){y — x) for all x,y. Define 



Vin,iy) = sup M ^{y)v{dy) - / ^(x)/i(dx) \ . (2.4) 

ip,5es U J J 

Weak duality gives that V{fj,, u) > P(/i, v) > T){p,, v). Note that there can be no uniqueness 
of the dual optimiser: if [tp, S) €z S then so is {ip + ax + b, 6 — a) . Hence we may choose any 
convenient normalisation such as V'(^o) = = 6{xo) for some xq £ M. 

In this article, in addition to requiring that // and ly are increasing in convex order we will 
make the following additional assumption. 

Dispersion Assumption 2.1. The marginal distributions /i and u are such that the support 
of rj := {fi — z^)+ is contained in an interval E and the support of j := (v — /i)+ is contained in 

Weak duality gives that V{^,v) > P(/i, i^) > D(/i,i/). Our goal in this article is to show 
that V{fJ,, v) = 'D(/U, v) and to give explict expressions for the optimisers p, ip and 6. 

Remark 2.2. Note that we make no other regularity assumptions on the measures /i and v. 
For example we do not require that /i and u have densities. In contrast, Henry-Labordere 
and Touzi [15j assume that /x has no atoms. This is also a simplifying assumption in part of 
Beiglbock and Juillet [3]. Conversely, the example in Beiglbock et al [2j which shows that the 
dual optimiser may not exist makes essential use of the fact that // has atoms. In our case. 



under Assumption 2.1 we show that a dual maximiser exists whether or not /i has atoms. 



Note that if// is less than or equal to v in convex order and Assumption 2.1 holds then E must 
be a finite interval. E may be closed, or open, or half-open. The rationale for Assumption 2.1 
will become apparent in the development of the results below. Let us, however, briefly point 
out that the assumption is natural in contexts most commonly encountered in mathematical 
finance. For instance, if the two distributions are increasing in convex order and log-normal, 



then Assumption 2.1 is trivially satisfied. 



Example 2.3. Let Z he a random variable which is symmetric about zero and which has density 
fz such that zfz{z) is unimodal on M-|_. For s < t let X = sZ and Y = tZ. Let fi and v be the 



laws of X and Y. Then Assumption 2.1 holds 



2.1 Heuristics and motivation for the structure of the solution 

Let B = {Bt)t>o be a standard Brownian motion and consider the problem of maximising or 
minimising E[|-B^|] over all stopping times r, subject to the constraint that E[-B^] = E[t] = 1. 
For the maximum, the solution is a two point distribution at ±1. This is consistent with the 
solution derived in Hobson and Neuberger |20j for the forward starting straddle, where in the 
atom- free case the solution to the problem of maximising E[|i^ — -'l^l] subject to X ~ //, y ~ z^ 
and E[5^|^] = X is characterised by a pair of increasing functions f,g with f(x) < x < g{x), 
such that, conditional on the initial value of the martingale being xq, the terminal value of the 
martingale lies in {f{xo),g{xo)}. Unfortunately, the condition that / and g are increasing is 
not sufficient to guarantee optimality and a further 'global consistency condition' ([20l p42]) is 
required. 

A solution to the problem of minimising E[|-Bt-|] over stopping times such that E[-B^] = 1 
does not exist. However, E[|-Bt-|] can be made small by placing some mass at in and a majority 
of the mast at 0. For instance, placing an atom of size ^^ at =bn and the remaining mass 
1 — l/2n^ at 0, we have E[|-Bt|] = l/n. The intuition which carries over into the minimisation 
problem for the forward starting straddle is therefore to move as little mass as possible. In 
other words, we expect that subject to the same constraints as above, K[\Y — X\] is minimised 
if y G {p{X),X, q{X)} for functions p, q with p{x) < x < q{x). 

2.2 Intuition for the dual problem: the Lagrangian formulation 

As in [20] we take a Lagrangian approach. The problem is to minimise J J \y — x\p{dx,dy), 
subject to the martingale and marginal conditions J^p{dx,dy) = v{dy), J p{dx,dy) = p{dx), 
and fy{x - y)p{dx, dy) = 0. 

Letting a(y), P{x) and 9{x) denote the multipliers for these constraints, define the La- 
grangian objective function L(x, y) = \y — x\— a{y) — f3{x) — 9{x){x — y). Then the Lagrangian 
formulation of the primal problem is to minimise over all measures p on M^ , 



p{dx,dy)L{x,y) + / a{y)v{dy) + / (3{x)v{dx). (2.5) 

For a finite optimum we require L(x, y) > and we expect that L{x, y) = for y £ {p(x),x, q{x)}, 



for a pair of functions (p, q), to be determined, see Figure 2.2 In particular, since L{x, x) = 0, 



we expect that /3 = —a. In this case we write L = L°'^ where 

L"'^(x, y) = \y-x\- Q(y) + a(x) - 9{x){x - y) (2.6) 

Assuming that a, 9, p and q are suitably differentiable we can derive expressions for a and 
9. The conditions L{x,p{x)) = 0, Ly{x,p{x)) = 0, L(x,q{x)) = and Ly{x,q{x)) = lead 




Figure 1: We have that a{x) + 0{x)\y — x| > a{y) with equahty for y £ {p{x),x,q{x)}. 



directly to the following equations for x £ E: 

X — p{x) + a{x) — a{p{x)) — 6{x){x — p{x)) = 

-l-a{pix)) + 6ix) = 

q{x) — X + a{x) — a{q{x)) — 6{x){x — q{x)) = 

l-a{qix)) + e{x) = 



(2.7) 

(2.8) 

(2.9) 

(2.10) 



Differentiating (2.7) and using (2.8), we obtain 



1 + a{x) - 9{x) - 0\x){x - p{x)) = 0. 



(2.11) 



Similarly, using (2.9) and (2.10), we obtain 



1 + a'{x) - e{x) - e'{x){x - q{x)) = 0. 



(2.12) 



Subtracting the second of these two equations from the first, it follows that 9'{x) = -r^x^j-x- 
Hence 9 is increasing and for some xq £ E, 



9{x) = 9{xo) + 



XO 



j{z) -p{z) 



dz. 



Now by adding (2.11) and (2.12) we have a'{x) = 9{x) + '^f-{2x-p{x) - q{x)). Then 



a{x) = x9{x) 



q{z)+p{z) 
XO <l{z)-p{z) 



dz + a(xo) — xq9[xq). 



Moreover, from the fact that q{x) is a local minimum in y of L(x, y) we expect < Lyy(x, q{x)) = 
—a"{q{x)). Hence a is concave at points y = q{x). 

Fix x' < x" . We must have that for y > x' , a{y) < a{x') + {y — x'){l + 9{x')), with equality 
at y = q{x'). We also know that for y > x" , a{y) < a{x") + {y — x"){l + 9{x")), with equality 
at 2/ = q{x"). Suppose that x" < q{x'). Then 



a(x') + (g(x') - x')(l + 9{x')) = a{q{x')) < a[x") + {q{x') - x"){l + 9{x")). 



Similarly, since x' < x" < q{x"), we have 

q(x") + {q{x") - x"){l + d{x")) = a{q{x")) < a{x') + {q{x") - x'){l + d{x')). 

Then adding and simplifying we find {q{x') — q{x")){9{x") — 9{x')) > 0. But 6 is increasing 
and thus q{x') > q{x"). Hence q is locally decreasing, in the sense that if x' < x" then either 
q{x') < x" or q{x') > q{x"). 

Now suppose we add Assumption |2.1[ The pair (x, q{x)) must be such that x & E, q[x) S E'^. 
Then, if x' < x" G £', we cannot have q{x') < x" . Thus q{x') > q{x") and q is necessarily globally 
decreasing. Similarly p is globally decreasing. 

3 Determining (^, q) 

3.1 A differential equation for p and q. 



Suppose that Assumption 2.1 is in force. Suppose that E has end-points {a, 6}, either of which 
may or may not belong to E. Then p : E >-^ (— oo,a] and q : E >-^ [6, — oo) are decreasing 
functions. 

Recall the definitions t/ := (/i — i/)"*" and 7 := (z^ — fJ,)^ . Our martingale transport of fi to i' 
involves leaving common mass {fi A i^) unchanged, and otherwise mapping r/ to 7. 

Suppose that t] and 7 have densities /^ and /^. Then we expect that p and g are strictly 
decreasing. Then {X ^Y,X < z) = {p{z) <Y < a)U {q{z) < Y). It follows that ioi z e E we 
have both 

/ frf{u)du= f^{u)du+ I f.,{u)du, (3.1) 

and, from the martingale property, 

pz pc<:i na 

/ ufjj{u)du= / ufj{u)du+ / ufj{u)du. (3-2) 

./a -^ q(z) -^ v{^) 



Then by differentiating equations (3.1 ) and (3.2 ) we can derive a coupled pair of differential 
equations for the pair (p, q) : 

, = Qix)-x frjjx) _ q{x)-x fp,{x)-fu{x) 

q{x) - p{x) f-yipix)) ^ q{x) - p{x) /^(p(x)) - fu{p{x)) ' 

, = x-p{x) fnjx) _ x-p{x) ff,ix)-fuix) 

q{x) - p{x) f-yiqix)) ^ q{x) - p{x) /^(g(x)) - fu{q{x)) ' 

where /^ and fi, are the densities of // and ly which we also assume to exist. Moreover we expect 
the boundary conditions p{a) = sup{x : Fj^{x) < F^(a)} and q{a) = sup{x : Fjj{x) < 1} (and 
similarly p{b) = inf{x : Fr^{x) > 0} and q{b) = inf{x : Frj{x) > F^(a)}). 

Example 3.1. Suppose fi ~ f^[— 1, 1] and v ~ C/[— 2,2]. We obtain p'{x) = fT_^A and 

^'(^) = q{x)-p{x) ■ ^^^ ^(^) = ^(^) -P{x) o.'^d J{x) = li^) +P{x)- Note that j'{x) = q'{x) + 
q'{—x) = j^/~\^ + ^~hfx) ~ ~^- ^y symmetry we expect that p{x) = —q{—x) and then j(0) = 
so that j{x) = —X. To derive h, note that q'{x)h{x) = p{x) — x and so {h'{x) + j'(x))/i(x) = 
j{x) — h{x) — 2x. Then h'{x)h{x) = j(x) — 2x = — 3x so that h{x)'^ = A — 3x^ for some constant 

A > 0. It follows that q{x) = {h{x) + j{x))/2 = \/q{0)'^ — |x^ — x/2. The boundary condition 
g(— 1) = 2 then gives q{0) = -v/3 and thus q{x) = — — ~|^ ~^ . Then also a second boundary 
q{l) = 1 is satisfied. Further, p{x) = -Vi2-3a:^- x_ 

6 



Example 3.2. Consider a slight modification of the above such that /i ~ C/[— 1,1] and v 
l\J\-2, -1] + |C/[1,4]. Then, forx G (-1, 1), p'(x) = -| (-f^^^) andq'ix) = -4 ( ^^iM. 



with boundary conditions p(— 1) = —1 and (?(— 1) = 4. This pair of coupled equations looks hard 
to solve. Nonetheless we will show using the methods of subsequent sections that 

p(x) = --(5 + ix + \/25 -8x- 8xA ; q{x) = -- ('s + 4x - 5\/25 - 8x - 8xA . (3.5) 

Except when it is possible to find special simplifications it appears difficult to solve the 
coupled differential equations for p and q. Indeed when the support of ly is unbounded, some 
care may be needed over the behaviour of p and q at the boundaries of E. More generally, if 
/x and ly have atoms then we cannot expect p and q to be differentiable, or even continuous. 
Instead we will take an alternative approach via the potentials of fi and z/, which leads to a 
unique characterisation of {p, q). 

4 Deriving the subhedge 

Motivated by the intuition of the previous section let Q be the set of pairs of decreasing functions 
(p : E >-^ (— oo, a], (7 : -E I— 7- [6, oo)). Using the monotonicity of p and q we can extend the domains 
of p and q to [a, b]. Note that we do not assume that p and q are differentiable or even continuous, 
or that they are strictly monotonic. Thus, in general p~^ and q~^ are set valued. Nonetheless 
we can define p~^ and qj^ to be the right-continuous and left-continuous inverses of p and q 
respectively. Then p~^ (y) = inf{x G E;p{x) < y} with the convention that p^ (y) = 5 if this 
set is empty. Similarly q^ (y) = sup{x G E; q{x) > y} and qj^ (y) = a if q{a) < y. 

Definition 4.1. Fix xq £ [a,b]. For {p,q) £ Q define 9 = Op^g : [a,b] — )■ M and a = ap^q : 

\a, 61 -^ M via 



eix) 



2dz 



q{z)-p{z)' 



I ^ f ^x-q{z)-p{z) f q{z)+p{z) 

a{x) = / — — — dz = xe{x) - / — —-dz. 

Jxo 1\z) - p{z) Jx^ q(z) - p(z) 

Extend these definitions to M by defining 6 = 6p^q : M — )■ R and ip = ipp^q : [a, 6] — )■ R via 

^{Vlli^)) X <a, 

d{x) = { eix) x£ [a,b], (4.1) 

6{p]^^{x)) X > b. 

(^{Pr^{^)) + {Pr^{x) - x){l - e{p]^^{x))) X <a, 
■0(x) = { a{x) xe[a,b], (4.2) 

Lemma 4.2. For x > b (respectively x < a) the value of tp{x) does not depend on the choice 
q2 £ q~^{x) (respectively pj^ (^p~^{x)). 

Proof. Let (^i)j=i,2 satisfy q{zi) = x. Then 

a{z2) + {Z2 - x){-l - 6{Z2)) - {a{z,) + (zi - x){-l - 6{zi))} = H ^^^^iM^z. 
But this expression is zero since q{z) = re on [zi, Z2\. D 



Example 3. 1 , Setting xq = 0, we have 

2 



1 



du 



Example 4.3. We calculate 9 and a in the uniform case for which we have found p and q in 
'ettim 

e{x) 

2 - V4 - a;2 



^du = —= I — : = —= sm 



and 



a{x) = x9{x) + 



ndu 



2x 



, ^ sin ^ f- I + 

Vl2 - 3u2 V3 V2/ 



^/3 



Recall the definition in (2.6) of the Lagrangian, now extended to functions ijj and 5 defined 



on 



L^'^{x, y) = \y - x\+ ip{x) - ip{y) + 6{x){y - x). 



(4.3) 



In Theorem 4.5 below we show that the Lagrangian based on (4.3) is non-negative in general 
and zero for x G {p{x), x, q{x)}. It follows that the pair {ip, 6) is a semi-static subhedge for the 
forward starting straddle, and hence that for any pair of random variables with E[l'|-'l^] = X 
and X ^ jj,, Y ^ v, 

E[\Y -X\]> I i^iyHdy) - [ i;{x)fi{dx). (4.4) 



Further, given a pair of strictly monotonic functions p : E >-^ (— oo,a] and q : E >-^ [6, oo), 
a probability measure ^ with density /^ and a sub-probability measure tj < fi with support 
contained in E and density /,,, we can construct a pair {X,Y = Yp^q^ri,p) via P(X G dx,Y G 
dy) = p*{x,y)dxdy where 



P*{x,y) 



q{x)—x 



fv\^)q(x)-p{x)hy=P(x)} y < ^' 
fii{x) - f-r,{x) y = X, 

. J'vi^) q{x)-p{x)hy=<i{^)} y>^- 



Certainly X ^ fi. Let z> = Vp^q^n.fj. ~ ^(X)- Then p*{x,y)dxdy G M.{p, v) and 
p*{x, y)\y - x\dydx = W.[\Y - X\] > P(/x, j)) > P(//, 9) > [ ijp^g{y)9{dy) - [ ijpjx) p{dx) 



and there is equality throughout, since p* is only positive when L'^P'1'°p-i = 0. In particular, 
V{p, v) = 'D{fi, jy) and there is no duality gap. 

Suppose p and rj are fixed, and consider varying p and q. Provided it is possible to choose 
p and q such that D = Vp^q^q^^ ~ v then we have 

V{p, ^)< I I P*i^^ y)\y - Adydx = E[|r -X\]= I ^p^g{y)u{dy) - f i;pjx)p{dx) < V{p, u). 

Weak duality gives V{p, y) > T>{p^ v), and hence the solutions to the primal and dual problems 
are the same. 

Remark 4.4. We have that for any (p, q) G Q 



inf E[\Y -X\]> / i^pJyMdy) - / i;pjx)p{dx) 

p&M{tM,u) J J 

In financial terms this means that any {p, q) £ Q can be used to generate a semi-static 
subhedge {ipp,q,Sp^q) and hence a lower bound on the price of the derivative. This bound is 
tight if, given Tg-call prices are consistent with Ftq ~ p, and Ti-call prices are consistent with 
Fti r^ u, then i/ ~ '^p,q,v,^J.■ 



In the case where Assumption 2.1 holds we show in the next section how to choose functions 
p and q such that Vp^q^rj^fi ~ i^ using a geometric argument. In fact, in cases where /i has atoms, 
genericahy there is no pair of functions (p, q) such that i'p^q^r],fi ~ i^ smd we have to work with 
multi-valued functions, or rather to re-parameterise the problem and use some independent ran- 
domisation. For the rest of this section we concentrate on showing that the choice of multipliers 



in Definition 4.1 leads to a non- negative Lagrangian. 
Theorem 4.5. For all x, y G M, 

L'f''^ = \y-x\+ ip{x) - ^{y) - 6{x){x - y) > 0, 
with equality for x £ E and y G {p{x), x, q{x)}. 
Proof. Suppose first that x G [a, ^]- If a^ < y < b, then 



q{z)-p{z) q{z)-p{z) 
-, / 2(,(z)-,) \ ^ 
\q{z)-p{z)) 

since q{z) >b>y. Now for a < y < x, 

L^'\x,y) = x-y + y(0(x)-%)) + £^|^|±^dz 

" dz (-1 + ^^ _ Pi^)+<li^) 

V q{z)-p{z) q{z)-p{z) 

^' , 2{y-piz)) 

dz—^ > 

q{z) -p{z) - 

since p{z) < a < y. 

Next, suppose y > b. Let j = q]^ (y)- Then "0(2/) = (y ~ t) + V'It) + ^(,l){y ~ l) and 
L {'y,y) = by construction. It follows that 

L'f''^{x,y) = {y-x)+i){x)-ij{-f)-{y-j) -6{j){y-j) + S{x){y-x) 
= a{x) - xe{x) - a(7) + 761(7) + ^-x + y{9{x) - 9{-/)) 



xo<i{z)-p{z) Jxoi{z)-pi^) Jx Jx q{z)-p{zj 

"^ dz 

-y- -—2{q{z)-y). 

q{z) -p(z) 

Now either q{x) < y or q{x) > y or q{x) = y. In the first case, since q is decreasing, x > 
1l (y) ~ T- Since for z G (7, x), q{z) < y we have L^' {x, y) > 0. On the other hand if q{x) > y 
then X < qj^ (y) = 7 and for z G (x,7), q{z) > y and again L^' {x,y) > 0. Finally if g(x) = y 
then X G q~^{y) and 9(2;) = y on (x,7) so that L"^' = 0. Similar arguments apply when y < a. 
Finally suppose x ^ [«,&]• We cover the case x < a, the case x > b being similar. For 
y <p]^^{x), 

< L^'-^C^J^^x), y) = V(P/jH^)) - V'(y) + Pe'{x) -y + 0[p-^\x)){y - PrH^))- 



But also L^'^{p^ {x),x) = 0, and thus 

il^ix) = i;{p]^^{x))+p^^{x) -x + e{pjl{x)){x-p'^{x)). 

Hence 

Q<L^^^{p-R\x),y) = x-y + ^{x)-^{y)-e{p-^\x)){x-y) 

< \x-y\+ V'(x) - V(y) - 9{p^\x)){x -y) = L^'^{x,y). 

For y > p]^ (x), 

L^''{x,y) - L^''{p-^\x),y) = 2{p-\x) - x) > 0, 

and hence L^'^{x, y) > 0. □ 

Corollary 4.6. L"^' {x,y) = for all pairs x G (a, 6) andy G [lini^j^^; ^l-z); hm^-f-^; g(z)]. Similarly 
L'^'\x,y) =Ofory£ [l[m^i^p{z),limz^^p{z)]. 

Proof. Suppose that q has a jump at x G -E and that y G [hm^j^^ q{z), liiaiz^x Qiz)]- If y is such 
that y / q{x) then we can modify q by defining q{x) = y and q{x) = q{x) otherwise. 

Write ip = tpp^q and 5 = 6p^q. Clearly, changing the definition of q at the single point x makes 
no difference to the definitions of tp or 6. Then, 

L^^Hx,y) = L'^Hx.y) = L^'\x,q{x)) = 0, 



the last equality following from Theorem 4.5 



D 

5 The general case 

The goal in this section is to show how to construct a martingale coupling, in the case where 



Assumption 2.1 holds. The aim is to construct appropriate generalised versions of p and q, 
allowing for atoms, in such a way that if X ~ /i, y G {p{X) , X , q{X)} and E[y|X] = X then 

Define D{x) = Cy{x) — C^{x). Then D is a non-negative function such that limu|_!.oo D{x) = 
and, because of Assumption |2.1[ D is locally convex on E'^ and locally concave on E. 

5.1 Case 1: /i A z/ = 

Suppose first that ^u A i^ = so that rj = fi and j = v are probability measures. 

Note that not both r] and 7 can have atoms at a. Let 7a = 7((— oo,a]) G (0,1). Then 
-fa = ^ax{D'{a-),D'{a+)}. 

Let F^^ : [0, 1] 1— ;■ [a, b] be the left continuous inverse of r/((— c«, x]) = r]{[a, x]). For u G (0, 1) 
define £u(x) = D{F:^^{u)) + (x — F^^{u)){'^a — u) — D{x) and 

£u{q) - D{p) 
(p(u) = sup . 

p<F-\u)<q l-P 

In regular cases, 4> is the slope of the line of which is tangent to both D on (00, a] and £u on 
[6,00). Let P{u),Q{u) be numbers such that P{u) < o, Q{u) > b and 

'^^"^ = Q{u) - Pin) ■ 
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P(u) 
y^D(P(u))++(u)(y-P(u)) 




f'iu) f'(v) QW Q(") 



Figure 2: The left figure shows, for a single value of u, how 0(it) is the slope of the line which 
is tangent to D{-) at P{u) and Sui') at Q{u). The right figure shows how as u increases <£'«(•) 
decreases and then P, Q and (j) also all decrease. 



It is clear from the construction, see also Figure [STU that </>, P and Q are decreasing functions. 
Define 

C[u)= mf 7a -n- (/)(«). (5.1) 

Then C is the slope of the line which is tangent to both Suix) = D(F^^{u)) + (x — F^^{u)){'ya — 
u) — D(x), defined over x < a and D(x) defined over x >b. 

Note that if P{u) is unique and D' is continuous at P{u) then D'(P(u)) = (t){u). Similarly, 
if Q{u) is unique and D' is continuous at Q{u) then E'^{Q{u)) = (j){u) which is equivalent to 
C{u) = D'{Q{u)). Since (f) and C, are decreasing functions and D'{P{u)) + D'{Q{u)) = 'ja — u, it 
follows that (j) is Lipschitz with unit Lipschitz constant. The following lemma shows that this 
holds in general. 

Lemma 5.1. (/) is a decreasing Lipschitz function and for < u < u < 1, (j){u) > 0('u) > 
(t){u) + (u- u). 

Proof. It is clear that £u{y) is decreasing in u. Then, for ii > u, 

, £u{Q{u)) - DiPju)) , SuiQju)) - D{P{u)) 

u) > ^rv—. ;rTT^ > :rTT^ ^tttt = 0{u) 



Q{u) - P{u) 
and hence (p is decreasing. Moreover, 
£u{Q{u)) - D{P{u)) 



Q{u) - P{u) 



(^) > 



Q{u) - P{u) 

^^ Q{u)-F-\u) D{F-\u)) - D{F-\u)) + {F-\u) - F-^){u){^a 



u) 



Q{u) - P{u) 



Q{u) - P{u) 



By construction the fraction in the middle term lies in (0, 1). Further, since D is concave 
on (a,b), D lies below its tangents and the numerator in the final term is positive. Hence 

(j){u) > 4>{u) + u — u. n 

Since (p is Lipschitz it is also absolutely continuous and we can write (p{u) = 7a + /g (j)'{w)dw 
for some function 6' . 
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Lemma 5.2. We have D{P{u))-P{u)(P{u) = f^ P{v)(t)'{v)dv, D{Q{u))-Q{u)C{u) = f^ Q{v)C'{v)dv 
and D{F^^{u)) - F-^{u){ja - u) = - £ F-^{v)dv. 

Proof. We prove the first statement, the other two being similar. Given {D(x);x < a}, for 
9 e (7a — 1, 7a) define the conjugate function Dc{9) = inf2;<a{^(a;) — 6x}. Then -Dc(O) = and 
Dc is concave so that Dc{6) = /q D'^{x)dx- Note that (j){u) is an element of the subdifferential 
of D at P{u) and hence 

D{Piu)) - Piu)^{u) = D,{<P{u)) = [ " D',{x)dx = - I (^'{v)D'M{v))dv. 

Jo Ju 

Finally note that D'^{(j){v)) = —{D')^^{(p{v)) = —P{v), almost everywhere. D 

Note that (j) has been constructed to solve 

4>{u){Q{u) - P{u)) = D{F-\u)) + (Q(n) - F-\u)){^a - u) - D{Q{u)) - D{P{u)). 

Using (j){u) + C,{u) = ^a — u we obtain 

D{P{u)) - P{u)4>{u) + D{Q{u)) - Q{u)au) = D{F-\u)) - F-\u){^a - u) (5.2) 

and hence 

P{v)(t)'{v)dv+ I Q{v)C{v)dv = - I F-^{v)dv. 

J u Ju 

Then we have both 0'(u) + C'{u) = —1 and P{u)(j)'{u) + Q{u)C,'{u) = —F~^{v) Lebesgue almost 
everywhere on (0, 1). Hence 

Q{u)-F'^{u) 

(j) (u) = — — ——— almost everywhere. 

Q{u) - P{u) 

Theorem 5.3. Let U and V be independent uniform random variables. Let X = F^^{U) and 
Y = P{U)LA,^^y, + Q{U)Iac^^^^ where 

_f ^ Q{u)-F-\u) \ 
"^(«-)-f - (Qiu) - P{u)) j ■ 

Then X has law rj and Y has law 7. 

Proof. It is immediate that X r^ rj. For y < a let P^ = infju G (0, 1) : P{u) < y}. Then 

/•I Q{u)-F-^u) 
JPrHv) Q(.u)-P{u) 

4>'{u)du = HPr^v)) = D'{y+) = 7((-oo,y]). 
A similar argument for y > h gives that y ~ 7. D 
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5.2 Case 2: /i A z/ 7^ 



If {i^i A i/)(]R) = 1 then we must have Y = X and ^\Y - X\ =0. Otherwise, let (/x A z^)(M) = 
K G (0, 1). Then r7(]R) = r]{E) = 1 — k = "){E^) = 7(M). Define probabihty measures fj, A u, r] 
and 7 by /I7w^(^) = (/U A u){A)/k, fj{A) = r]{A)/{l - k) and j{A) = j{A)/{l - k). 
The following result follows immediately from Theorem |5.3[ 



Theorem 5.4. Let U and V be independent uniform random variables, and let Z be a Bernoulli 
random variable with parameter k. Conditional on Z = 1, then set Y = X ^ fj, Av. Otherwise, 
conditional Z = 0, let X = F^ (U) and Y = P{U)lA,^y. + Q{U)Ij^c where 



A 



{u,v) 



V < 



Q{u) - Fz±{u) 



[Q{u)-P[u)) 



where P and Q are defined as in Section [577] for the disjoint measures rj and 7. 
Then X has law /x and Y has law v. 

5.3 Strong duality 



Theorem 5.5. Suppose fi < u in convex order and Assumption 2.1 holds. Then V{fi,i' 



P(//, z^) = P(//, z/). Moreover, the optimal coupling is given by the construction in Theorem 



5.. 



and the optimal semi-static hedging strategy is given by the expressions in (4-1) and (4-2), for 
{p{x),q{x)) = {P{F^{x)),Q{F^{x))). 



Proof. By Theorems |5.3| and |5.4[ we have a candidate coupling. Moreover, by Theorem 4.5 
and Corollary 4.6, it follows that L'^' = wherever the candidate coupling places mass, and 
hence, for this choice of joint law for {X,Y) and (^,(5) = {ipp,q,6p,q) we have K[\Y — X\] = 
f ^{y)v{dy) - f i/j{x)^{dx). Then P(^, u) <E[\Y- X\] = f ^{y)y{dy) - f ^{x)fi{dx) < V{fi, v) 
and by weak duality there is equality throughout. D 



6 Examples 

6.1 The case with densities 

If // and V have densities then it is not necessary to introduce P and Q and instead we can use 



the functions p{x) = P{Frj{x)) and q{x) = Q{Fr^{x)). Then from (5.1) and (5.2) we find that 
the pair (p, q) are characterised by 



D'{x) 
D(x)-xD'(x) 



D'{q{x)) + D'{p{x)), 

D{q{x)) - q{x)D'{q{x)) + D{p{x)) - p{x)D'{p{x)), 



(6.1) 
(6.2) 



for X £ E. 



Example 6.1. Recall Example 3.1 and suppose fi ~ U[—l, 1] and v ~ C/[— 2,2]. 

D(x) = 



Ix^ + Ix + l 



\J2 



1^2 1 1 
8 2 



-2< X< -1, 
-1 < X < 1, 



X- - ^x^\ \<x<2. 



By (g. 1\ ), for —1 < X < 1 we have —x = q{x) + p{x). Then, by (6.2) for x G [—1, 1] we have 

D{x) - xD'{x) = D{q{x)) - q{x)D'{q{x)) + D{-x - q{x)) + (x + q{x))D'{-x - q{x)). 
After some simplification we find q{x)'^ + q{x)x + x^ — 3 = 0, whence q{x) = ^'^^ ^ ~ ; ^'^'^ 
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Example 6.2. Recall Example 3.2 and suppose that fj, ~ L'^[— 1,1] and u ~ |[7[— 2,— 1] + 
I [/[1, 4]. Then, 



D{x) 
with D = outside [—2,4]. 



4(^ + 2)2 

16 "•" 8 i^ 



-2<x < -1; 
-Kx < 1; 
1 < a; < 4, 



Then (6.1) and (6.2) yield 

4x + 5p + g + 5 = 0, 



4x^ + bp^ + q^ = 25. 



Eliminating q and solving for p gives the first expression in (3.5), from which the second follows. 

6.2 The general case with atoms 

The full power of the analysis is illustrated in the case where /i has atoms or u has intervals 
with no mass. Then there is no hope of constructing the optimal martingales via an approach 
using differential equations. 

Example 6.3. Let ^ he any measure supported on E = [—1, 1] and suppose v has support E'^ 
and density v{dx) = \x\^'^dx for \x\ > 1. Then for \y\ > 1, D{y) = i^r-r, D'{y) = — 2~n" ^'^^ 

D{y) - yD'{y) = ^. FuHher, Ja = h 

Since ^ is not yet specified we do not have an expression for D on E. Nonetheless, for 
9 E (—5, 5) we can define the conjugate Dc{0) = sup_i^^^i{D{x) — x9), so that Dc{\ — u) = 
D{F^\u))-E~^{u){\-u). 

Then P = P{u) and Q = Qiu) solve 



^{u) = D'{P)=8'^{Q) 



£u{Q) - D{P) 
Q-P 



and hence 



1 



1 



2P2 \^2 

Then, from the first equality 



u + 



Dr 



u) + 



u)Q - D{Q) - D{P) 



2Q2 



Q-P 



P2 Q2 



l-2u 



(6.3) 



and 



^'2-^+20^ 



P 



2P2 



'2 ) \1 ) 2Q 2P' 



which simplifies to 



1 



\F\^^ = '''\2 



u 



(6.4) 



It is straightforward to solve (6.3) and (6.4) for 1/P and 1/Q and we find 

2D,{l-u) 2Dc{l-u) 



P{u) 



D,a-u)^ + 2a-u) 



Q{u) 



D,a - u)^- 2a -u) 



For example, suppose ^jl ~ 5'^-i/2 + |f^[Ojl]- Then, for x G ( — 1,— g] D{x) = 1 + x/2; for 
x G (-i 0), D{x) = 3/4; and for x G [0, 1], D{x) = (3 - a;2)/4. Further, for < u < 1/2, 



Piu) 



(2-u) 



2 - 3n + n2/4 
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Qiu) 



u) 



u + u'^/i' 



whereas, for 1/2 < u < 1, 



P(n) 



2{1 -u + u"^ 



2 - 4tt + 3w2 - 2u3 + u^ 



Qiu) 



2{1 -u + u^) 
3n2 - 2n3 + u^ ' 



Example 6.4. Suppose fi ~ U{—1,1) and v ~ C/{— 2,— 1,3}. T/ien D is zero outside (—2,3) 
and D{x) = '^ for -2 < x < -I, D{x) = \ + ^ - ^ for -l< x <l and D{x) = 1 - f for 
1 < X < 3. 

It is clear that q{x) = 3 and p{x) = — 1 for x < x* and p{x) = —2 for x > x* , where 
X* £ (—1,1) is to be found. 



Then, see Figure 6.4 • x* is the solution to 

{D{x) + (y - x)D'{x) - D{y)}\^^^ = {D(-1) + (y + I) D' {-!-)} \^^.^ 

and we find x* = 3 — 4 



3' 




Figure 3: The point x* is such that the tangent to D{-) at x* meets the hne y = (x + 2)/3 (ie 
the Hnear extension of y 1— )• D(y) taken over (—2, —1)) at x-coordinate 3. For x > x* the hne 
joining (— 2,L'(— 2) = 0) to (3,iS^(^.)(3)) hes below D on (— 00,— 1) and above £f{x) on (l,oo). 
Hence, for (x > x*) we have (p(x) = — 2,g(x) = 3). For x < x* (not shown) the hne joining 
{—1,D{—1) = 1/3) to (3, £^^(a.)(3)) hes below D on (—00, —1) and above £f{x) on (1, 00). Hence, 
for (x > X*) we have (p(x) = —1, ^(x) = 3). 



7 Concluding remarks 

7.1 Upper bound 

We show below that in the special case where i) (/i A i^)(M) = 0, ii) /x has support on an interval 
E and iii) 1^ has support in E'^, a modification of the above methods yields the martingale 
coupling which maximises K[\Y — X\]. Moreover, the optimal transport plan and the super- 
hedging strategies can also be given explicitly as for the lower bound case. Recall that Hobson 
and Neuberger [2D] give an existence proof, and show that there is no duality gap, but do not 
give a general method for deriving explicit forms for the optimal coupling. 

If // is a point mass at x, then the problem is trivial and E[|^ ~ x\] = J \y — x\i'{dy). 
Otherwise, proceeding as in Section 2.2 but noting that now we expect L < 0, /3(x) > —a{x) 
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and, conditional on X = x, mass at two points only, we can find equations for nionotonic 
increasing functions (now labelled g and h, instead of p and q). Then a is convex at points 
y = g{x) < X and y = h{x) > x. We find also that g and h are increasing. This general 
structure is evident in [20] . 

Now we add the following assumption: 

Strengthened Dispersion Assumption 7.1. The marginal distributions /x and v are such 
that {fj, A i^)(K) = 0, the support of fi is contained in an interval E and the support of i' is 
contained in E^. 

As before, suppose E has endpoints {a, b}. Then necessarily g : E >-^ {—oo, a] and h : E h^ 



[b, oo). In the case with densities we can write down equations similar to ( |3.1[ ) and (3.2): 

fg{z) fh{z) 



f^{u)du = I fu{u)du+ I fu{u)du (7.1) 

7 — oo J b 

r-g{z) r-h{z) 

uf^{u)du = / ufy{u)du+ I ufy{u)du. (7-2) 

la J — oo Jb 

The defining equations for g and h become C'^{x) = Cl{f{x)) + Cl{g{x)) — Cl{b) and 

C^ix) - xC'^ix) = a(/(x)) - f{x)CUf{x)) + C,{g{x)) - g{x)Cl{g{x)) - [C,{b) - bC'M]. 

This time the equations do not simplify into expressions involving D{x) = Cu{x) — C\i{x) alone. 
By analogy with the lower bound case, and now allowing for measures which are not abso- 
lutely continuous, introduce for u G (0, 1) and y > b, 



Fu{y) = C^{F-\u)) + {y- F-\u)){u - 1) + C,{b) + {y - b)CUb) - C,{y) 



and define 



ip{u) = sup 

g<a<b<h "- " 5 

Then the suprema is attained at G{u) < a < b < H{u) and 

f{u) = T7r-\ — 7v—\ = C^^C^M) = ^u{H{u)) 

H[u) — G(u) 

assuming these last two derivatives exist. 



Theorem 7.2. Suppose Assumption \ 7. 1\ is in force. 

Let U and V be independent uniform random variables. Let X = F~^{U) and Y 

G^^)^A,,^,,+H{U)I^c^^^^ where 

~ \ Biu)-F-\u) \ 

^"'^^ " 1 - (^M - G(n)) / • 

Then X has law fi and Y has law v. 

Moreover ^[\Y - A|] = suppg_^(^ ^,) S S \y - Ap{dx, dy). 



Proof. The proof of the first part of the theorem follows the proof of Theorem 5.3 The fact 



that this coupling yields the maximal value of J/ |y — x\p{dx,dy) follows from analogues of 



Theorem 4.5, Corollary 4.6 and Theorem 5.5 D 
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Example 7.3. Suppose fi r^ U[-l,l] and i/ ~ ^t/[-3, -1] + ^C/[l, 3]. T/ien C^(x) = (1 -x)V4 
for X E ( — 1,1) and C^{x) = x^ = (— x) V otherwise. Further, Cy{x) = ^^ „ ' x for 

-2,<x< -I, Cy{x) = 1 - f /or -1 < X < 1, Cy{x) = ^^^ /or 1 < x < 3 and C^{x) = x" 
otherwise. 

In this case g{x) = G{F^{x)) and h{x) = H{Ffj_{x)) are well defined. Define J-x{y) = 
Tp (x){y)- Then J-x{y) = g (-'^ ~ ^■^'^ ~ 2y + 4xy — y^) and it is easy to check that the pair 
{g{x) = X — 2, h{x) = X + 2) solve 

Hence suppg_/\^(„ u) 1 1 \y~x\pidx, dy) is attained by the coupling Y = X + Q where Q is uniform 
on {±2} and is independent of X . 

For this example there is a simple proof of optimality. Jensen's inequality gives ¥i[{Y—X)'^] > 
(E[|5^ — ^|])^j CLnd the left hand side is independent of the martingale coupling and equal to 
j y'^i>{dy) — f x'^fi{dx). Hence for our example, K[\Y — X\] < 2 for any martingale coupling, 
and K[\Y — X\] is easily seen to be equal to 2 for the coupling of the previous paragraph. 

Remark 7.4. Unlike in the lower bound case, it is not the case that if (// A i^)(M) > then 
the optimal construction is to embed the common part of the distribution by setting Y = X 



and then to embed the remaining part by the style of construction given in Theorem 7.2 This 
explains why we need the strengthened dispersion assumption. 

7.2 The situation when the Dispersion Assumption does not hold 



Return to the lower bound and consider the case where Assumption 2.1 fails. For simplicity 
we assume that // and v have densities. Then we expect there to be functions p{x) < x < q{x) 
such that conditional on X = x,Y £ {p{x),x,q{x)}. 



From the analysis of Section 3.1 we expect that if x' < x" then p{x") ^ {p{x'),x') and 
q{x') ^ {x" ,q{x")). This condition is necessary for optimality, but not sufficient and there may 
be several candidate martingale couplings characterised by pairs {p, q) with this property. A 
further condition is necessary, akin to the 'global consistency condition' of |2Uj . 

Example 7.5. Suppose fj, ~ U[—l, 1] and v ~ ^C/[— 2, —1] + ^i^o + jf^il; 2] where 5q denotes a 
point mass at zero. Suppose the mass on [—1,0] is mapped to [c, —1] U {0} U [d, 2] where —2 < 
c < — 1 and 1 < d < 2. Then mass and mean considerations imply ^ = — 4—^ + + ^ ^ ' and 
-y\ — — 4-^ 2 + 4 2 where 9 is the mass mapped from [—1,0] to zero, so < 9 < ^. 
These equations simplify to c + d = 49 — 1 and (P' + d'^ = 7. Then d = - — ^^ + — 2 ^^'^ 

c = - — 2"^ — 2 ■ ^^ particular there is a parametric family of martingale couplings for 

which p and q have the appropriate monotonicity properties. 

In this example we expect by symmetry that the optimal solution has 9 = ^, and this can 
indeed be shown to be the case. In general further arguments are required. 

7.3 Multi-timestep version 

The analysis of the lower bound extends in a straightforward fashion to multiple time-points. 



provided an analogue of Assumption |2 . 1 1 holds . 

Let To = < Ti < ... T„ be an increasing sequence of times and suppose the marginal 
distributions /ij of a martingale M are known at each time Tj. Consider the problem of giving 
a lower bound on 5 = E[Ej=i \Mt, - Mt^_^\]. 
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Figure 4: A viable pair (p, q), in the relatively simple case where the set with fi — v > consists 
of two intervals. 

For the problem to be well-posed we must have that the family (/ij) is increasing in convex 
order. Let r/j = {fj-i-i — m)^ and 7i = (fj-i — ^j_i)+. Suppose further that r/j has support 
contained in an interval Ei and ji has support contained in Ef. Then the methods of the paper 
apply, and we get a tight lower bound on S of the form S > X^"=i V{fii-i, fii). 



7.4 Model independent bounds and mass transport 

There have recently been several papers, including Beiglbock et al {2|, Beiglbock and Juillet 
[3] and Henry-Labordere and Touzi [T5] which have made connections between the problem of 
finding optimal model-independent hedges and the Monge-Kantorovich mass transport problem, 
and between the support of the extremal model and Brenier's principle. 

Beiglbock and Juillet [3] introduce a left-curtain martingale transport plan. This coupling 
has many desirable features, and using the methods of Henry-Labordere and Touzi [15j it is 
possible to calculate the form of the coupling in several examples. The motivations of Beiglbock 
and Juillet [3] in introducing the left-curtain coupling are at least threefold: firstly it is rela- 
tively tractable; secondly they use it to develop an analogue of Brenier's principle; and thirdly 
this coupling minimises E[c(X, y)] for any cost functional c{x,y) = h{y — x) where the third 
derivative of h is positive. (Note that this manifestly excludes the case h{z) = \z\ which is the 
topic of our study.) This optimality result was extended by Henry-Labordere and Touzi [15] to 
include all functions c which satisfy the Spence-Mirrlees type condition Cxyy > 0. (One example 
which is natural in finance is the payoff c{x, y) = —(log y — log x)'^ which arises in variance swap 
payoffs, see also [18J.) 

In addition to the clear financial motivation in terms of forward starting straddles, the payoff 
|y — x| is a natural object of study since it is the original Monge cost function in the classical 
(non-martingale) setting, see Riischendorf [23]. In that setting the non-differentiability makes 
this payoff a relatively difficult functional to study. 

The approach in [3] is based on a notion of cyclic monotonicity in a martingale setting. 
This has the advantage of providing existence results in a general setting but is typically not 
amenable to explicit solutions. In contrast, the Lagrangian approach employed in this article 
leads to tractable characterisations of the optimal coupling. A further contribution of this study 
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is that unlike both [3] and |15j who assume that fj, (at least) is atom-free, we indicate how to 
deal with atoms in /i. We find a geometric representation for the optimal coupling, which can 
then be applied to arbitrary distributions. It is likely that these ideas can be applied more 
generally. 
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